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Abstract 

The "correlated-projection technique" has been successfully applied to 
derive a large class of highly non Markovian dynamics, the so called non 
Markovian generalized Lindblad type equations or Lindblad rate equa- 
tions. In this article, general unravellings are presented for these equa- 
tions, described in terms of jump-diffusion stochastic differential equations 
for wave functions. We show also that the proposed unravelling can be 
interpreted in terms of measurements continuous in time, but with some 
conceptual restrictions. The main point in the measurement interpreta- 
tion is that the structure itself of the underlying mathematical theory 
poses restrictions on what can be considered as observable and what is 
not; such restrictions can be seen as the effect of some kind of super- 
selection rule. Finally, we develop a concrete example and we discuss 
possible effects on the heterodyne spectrum of a two-level system due to 
a structured thermal-like bath with memory. 

PACS: 42.50.Lc, 03.65.Ta, 02.50.Ey 

Keywords: Lindblad rate equation; unravelling; non-Markovian stochastic 
Schrodinger equation; heterodyne detection; direct detection; quantum 
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1 Introduction 

Open quantum system theory concentrates on the study of the time evolution 
a quantum system in contact with an environment; in particular, this theory 
aims to describe phenomena such as decoherence, relaxation, emission of light, 
evolution of entanglement [TH5]- Starting from the Hamiltonian approach de- 
scribing the coupled evolution of the quantum system and the environment, the 
reduced evolution of the quantum system is obtained by tracing out the degrees 
of freedom of the environment. This allows to describe the time evolution of 
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the open system in terms of its density matrix ps (t) with the help of a quantum 
master equation. Invoking standard assumptions as weak coupling limit and 
Born-Markov approximation, one can derive the Markovian quantum master 
equation [3HS], with infinitesimal generator in the Lindblad form [JJJE]. This 
approach, called the Markovian approach, is physically based upon the absence 
of memory effects in the action of the environment. This is a good and useful 
assumption in several physical examples, namely in quantum optics [IH5]- 

However, such assumptions are not valid in general and in many physically 
important cases the description of a reduced quantum evolution requires a non- 
Markovian approach involving strong and long memory effects. For example, 
situations with strong coupled systems, entanglement and correlation in the ini- 
tial state, finite reservoirs. . . need to be described by non-Markovian dynamics. 
Different techniques, such as the Nakajima-Zwanzig projection technique, the 
time-convolutionlcss operator technique, random Lindblad operator, random 
functional equations have been developed to derive non-Markovian quantum 
master equations [3ll4llMl4]. Recently, the concept of correlated projection 
technique has been used in order to describe a non-Markovian generalization of 
Lindblad type master equations (or Lindblad rate equations) [T5HT7] . This ap- 
proach has been successfully applied to describe non-Markovian models: struc- 
tured reservoirs, two-state systems coupled with energy bands [TT)H2"5] . . . 

An active line of research concentrates on the study of the behaviour of 
the solutions of these equations (thermalization, return to equilibrium, deco- 
herence,...). But even in the Markovian case, the quantum master equations 
remain often of a formal interest. In particular, most of the equations cannot 
be solved analytically and involve a large number of parameters which prevent 
numerical simulations. Concerning the numerical aspect, a powerful approach 
is the theory of "stochastic wave function unravelling". This consists in con- 
structing a stochastic differential equation for a wave function ip(t) such that 
M[\tlj(t))(tp(t)\} = ps{t). Then, by taking the average of a large number of re- 
alizations of ip(t) one reproduces the solution of the master equation. This 
has been applied in many Markovian situations [2,4,26,27]. Concerning the 
non-Markovian framework, different extensions of this approach have been de- 
veloped [2"8"H30] (there is no general and common approach). 

In the Markovian case the stochastic unravelling of the master equation has 
not only a technical usefulness, but it can be also interpreted in terms of mea- 
surements in continuous time; often the name of quantum trajectory theory is 
used [2,27 . In particular, for quantum optical systems the stochastic formula- 
tion is used to describe direct, heterodyne and homodyne detection. However, 
in the non-Markovian setup the notion of quantum trajectories as well as the 
measurement interpretation are still highly debated 8.9. 28 . 30T433] . 

For the non Markovian generalization of Lindblad type master equations 
[T5HTT] , only particular unravellings have been presented . In this article, 

we aim to present a general approach to obtain unravellings for this type of 
equations and to show that in this case an interpretation of the unravelling in 
terms of measurements in continuous time is possible. Our approach is based 
upon the general technique used to unravel Markovian Lindblad equations. In 
particular, our results include and generalize the previous results |34j . However, 
we have an important conceptual difference from the Markovian case. We are 
assuming that the structure of the bath responsible of the non Markovian be- 
haviour is not observable and this makes unobservable some of the components 
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of the noises introduced in the unravelling. 

The article is structured as follows. In Section [3J we describe the Lindblad 
rate equations. In Section^ we present the jump-diffusion unravellings of these 
equations. In particular, we derive non Markovian generalizations of stochastic 
Schrodinger equations. The stochastic master equations and the measurement 
interpretation are given in Section 2J In Section [SJ we construct a concrete non 
Markovian model (a two level system in contact with a structured environment), 
we present a possible unravelling, and we show possible effects of the non Markov 
dynamics on the heterodyne spectrum. 

2 Non-Markovian generalized Lindblad-type 
master equations 

In this section we introduce the non Markovian Lindblad-type master equation 
which we are interested in. These equations can be obtained by the application 
of the correlated projection technique and are sometimes called Lindblad rate 
equations [Tr)llT5H2"Tll2"3T - l2"5] . For any separable complex Hilbert space Oi we 
denote by L ( J£) the space of the linear bounded operators on J£, by T( J£) the 
space of trace class operators and by S(!K) the set of statistical operators (a 
statistical operator is a trace class, positive operator with trace 1). 

Let !Ks denote the Hilbert space representing the open system. The gener- 
alized master equation we consider is the evolution equation 

i*(t) = -i[H\ m {t)] + E E (K k Vk(t)K k * \ {RfR^Mt)}) (1) 

for the vector (j]i(t), . . . , ?7n(t)) with components in T(!Kg). The quantities 
H l , R 1 ^ are system operators which we take to be bounded for mathematical 
simplicity and A is a finite set of indices. 

Assumption 1. IK5 is a complex separable Hilbert space, H l = H l * e £(!Ks), 
Ra S £>CKs), k,i = 1, ... ,n, a € A. The initial condition of Eq. (fTJ) has the 
properties 

n 

Vi(0) 6 T(JC S ), 77,(0) >0, ^Tr Ms {r ?i (0)} = l. (2) 

i=i 

Remark 1 . Equation (fTJ) preserves the properties © at all times [TB] ; then, we 
interprete as system state the statistical operator 

n 

Vs(t)=^2vi(t)- (3) 

i=l 

The proof of the positivity preservation property of Eq. (JTJ) is very instructive 
and goes through the embedding of the dynamics {?7i(0)} 1— > {rji(t)} into an usual 
Lindblad dynamics in an extended state space [16j . Let us consider the enlarged 
space 'K = !Ks (8> C™. Let {e.^, i = 1, . . . ,n} be a reference orthonormal basis of 
C n . Let us introduce the block diagonal operator r)(t) on 5f by 

n 

i=i 
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and set 

n 

H = ]T H* ® \e t ) (a | , S« = ® |ei> (e, | . (4) 
»=i 

Then, we get immediately from ([1]) the evolution equation for fj(t): 

^ t fj(t) = cm)] = -mm] + EE - \ {sTs%,m}) ■ 

a€Ai,j = l ^ ' 

(5) 

For block diagonal initial conditions the two equations ([1]) and ([5]) are completely 
equivalent. Let us note that the linear map C is explicitly in the Lindblad form, 
so that the maps fj(Q) i— > fj(t) and {?7i(0)}j™ 1 i— > {f?i(i)}j=i are completely positive 

( CP )' 

In spite of the construction above, the index i is not interpreted as a quan- 
tum degree of freedom, but as the value of a classical observable. In typical 
applications the index i labels the energy bands of a structured environment 
[T51H6, 21,22,34;. A vector of operators with the properties © can be seen as a 
classical/quantum state. If we set Pi(t) = Trjt s {?7i(i)} and fji(t) = i]i{t)/pi(t), 
we have fji(t) € S(Jts), p%(t) > 0, X)ILift(^) = 1- I n quantum information the 
set of probabilities and statistical operators {pi(t), f]i(t); i = l,...,n} is called 
an ensemble and it is completely equivalent to the vector (j]i(t), . . . , rj n (t)) 
[3S1[37]. In this setup the system state is known as average state and it 
does not contain the information on the classical label i. Equation ([T]) gives a 
memoryless evolution for the ensemble {pi(i), fji(t); i = 1, . . . , n}; it is the evo- 
lution of the system state rjs(t) which is non Markovian. 

In the Markov case it is well known how to construct general unravellings of 
a master equation and how to give a measurement interpretation to them. So, to 
have an usual master equation in Lindblad form extending our non Markovian 
dynamics ((TJ) is a good starting point for the whole construction. However, Eq. 
(JS|) is not the unique extension of ([T]) and here we give another extension which 
is in some sense more convenient as starting point. The possible extensions 
depend on having or not the condition Kg oc 5^-; so, we put in evidence some 
diagonal terms. 

Assumption 2. Let us take A = {—mi, . . . , —1, 1, . . . ,1712}, and R l l a = 5ijL l a , 
for a = 1, . . . , mi. 

With this assumption Eq. ((TJ becomes 

-nViit) = K.i(r]i(t),...,T] n (t)), (6a) 



mi / 1 \ 

£i(ri, • • • , T n ) := -i[H\ Ti\ + E Li « T * L « - o { L «* L «' T *} 

+ E E (K k T k Rf - \ {RfRt\n} . (6b) 

a=l k=l \ ~ / 

By using the operators (0]), we define the new operators 

n n 

^ = E^» E E i « 8 l ei ^l> a=l mi, (7a) 

i,j=l i—1 
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n n 



S^ = J2 = £ R l i ® 1*1, P = 1, . . .,m 2 , 



(7b) 



i=l i=l 



and the Lindblad map C: Vr e T(5-{), 




(8) 



a=l 7=1 ^ ' 



Then, we consider the Markovian quantum master equation 



(9) 



with the initial condition 



?7(0) e §(JC), 



Tr c , {77(0) (l®|e i )(e i |)} = » 7i (0). 



(10) 



Remark 2. Let us use a subscript i to denote the i-th block on the diagonal 
of any trace-class operator, i.e. n = Trc« {r (1 (g) |e^) (e^ |)}- It is easy to check 
that 



and, so, both the master equations ([5]) and ^ reduce to the same Lindblad 
rate equation © for the blocks on the diagonal, while they are different for the 
off-diagonal blocks. Being equal at time t = due to (JTUJ), we have that the 
blocks on the diagonal of rj(t) are exactly the quantities rji(t) satisfying Eq. ([6]). 

Another way to describe the situation is to say that there is a superselec- 
tion rule and only block-diagonal observables are permitted. Then, statistical 
operators with the same blocks on the diagonal are equivalent and represent 
the same physical state. In this sense the two master equations ((5J) and © are 
physically equivalent. 

It is worthwhile to note that the operator C can always be written in the 
form C. It is enough to change the meaning of the subscript in the operators 
R l J or L l a in such a way that it includes also the index i. Then, given two triples 
(i,j,a) and (i',j',a'), we have that i 7^ i' => a 7^ a' (the same holds for two 
couples (i,a) and (i',a')). In this way, in the sums in Eqs. Q only one term 
survives and C = C. So, there is no loss of generality in considering only the 
master equation Q; the other case is always included, eventually at the price 
of a renaming and reordering of the indices. 

It is useful to formalize the framework we have presented in terms of normal 
states on W-^-algebras and of CP dynamics. 

Remark 3. Let G(X; L(J~Cs)) be the W^-algebra of the functions from X = 
{1,2,..., n} into £(^Ks) [3S1[35]. By natural identifications we have C(X; C) ~ 
C n and e(X;L(M s )) ~ £(Ks) ® C n , so that a £ e(X; £(5C S )) means a = 
(di,...,a n ), cij e £(?C); then, ||a| = max^gx ||%||- The predual space of 
e(X;£(5C s )) is S(X;T(JC S )) ~ T(5C S ) ® C n , so that r G e(X;7(K s )) means 



t = (ti,...,t„), rj G T(5C); then, Hr^ = £\ eX Htj^ = E"=i Tr M s (v 77 "/ 7 ^"}- 
In a natural way a and r can be considered as block-diagonal elements of £(JC) 
and T(!K), respectively: a ~ Sj=i a j ® l e j)( e jL r — Ej=i r i ® l e i)( e jl- 



£[r]j = £[r]j = ^(n,. . . ,r„), 



(11) 
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Remark 4. Equations © and © define a CP quantum dynamical semigroup 
T(t) on T(JC). Then, we define the projection V : T(5C) -> C(X; T(5C S )) C T(5C) 
by ("P[t]).,- = Trc"{r(ll ® l e i) ( e j 1} • The dynamics associated to the Lindblad 
rate equation © turns out to be V o TTi) I _/„ \ ; it is CP and Markovian. 

Finally, we define the projection V s ■ e(X;T(M s )) ->■ T(J£) by P s [r] = £V r,. 
The CP dynamics giving the system state (131) is Vs ° "P ° T(f)L/„ a and 

1 6 ^X;T(:Hs)J 

it is this dynamics which is non Markovian. 

3 Unravelling of non Markovian Lindblad-type 
master equations 

In this section, we derive a general form of jump-diffusion stochastic differential 
equations (SDEs) for wave functions in the enlarged space "K = "K s S3 C™ which 
provide unravellings of the Lindblad rate equations ©■ Having at hand the 
usual Markovian master equation ©, we adopt the usual approach pTIIMIHO] 
of stochastic Schrodinger equations in the Markovian case. This method is based 
on classical stochastic calculus (see for instance Refs. [4T[H2] and [3JJ Appendix 
A]) and the notion of a posteriori states [43ti46] . 

The key point of the theory is the construction of a linear and a non-linear 
stochastic Schrodinger equation (SSE), connected by a normalization and a Gir- 
sanov transformation, and, then, of the linear and non-linear stochastic master 
equations. The non-linear SSE is the key starting point for numerical simula- 
tions of the solution of a master equation, while the possibility of passing to 
linear equations is fundamental for the possibility of giving a measurement in- 
terpretation to the whole construction without violating the rules of quantum 
mechanics. Finally, the non-linear stochastic master equation gives the a pos- 
teriori states, the conditional state to be attributed at the system at time t, 
knowing the results of the measurement up to time t. 

3.1 The linear stochastic Schrodinger equation 

We consider a filtered probability space (f2,3 r , (3t),Q), satisfying the usual hy- 
potheses [271 Appendix A]. On this space, we consider d\ + c?2 x n independent 
standard Wiener processes W a , Wl (a = 1, . . . , d\ < mi; j3 — I, . . . , d 2 < m 2 ; 
j = 1, . . . , n) and (mi — d±) + (m 2 — cfe) x n independent standard Poisson point 
processes N a of intensity X a > and Ni of intensity A \ > (a = d\ + 1, . . . , mi ; 
(3 = di + 1, . . . , m 2 ; j = 1, . . . , n), also independent of the Wiener processes. 
All these processes are adapted and Wj(i), N a (t) — X a t and N 3 a {t) — X J a t are 
(3 r t)-martingales, under the reference probability Q [HJ25]. The trajectories 
of the Wiener processes are taken to be continuous and the trajectories of the 
Poisson processes continuous from the right. We set also 

mi m2 n 

E A «+ E E AJ - 

a— a— c?2 + l j — 1 

Now, on (5Ft),Q), we consider the following SDE for an J{-valued 
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process: 



dC(t) = [K + - C(*-)d* + E V a ((t-)dW a (t) 



a=l 



d 2 n 



££^C(*-)dW*(t)+ E (-^=^-l)C(t-)diV a (t) 



Q— 1 fc— 1 a— + 1 

m2 n 



C( = d2+1 fe=l 

where the operator in the drift part is given by 

mi mr> n n 

* = - 2 E v °<* v °< - 2 E E = E^' ® i e ^i< 

a— 1 a— 1 fe=l j=l 

1 mi rri2 ii 

Kj = 2 E L * L i - 2 E E ■ 

a— 1 a— 1 A:— 1 

By using the decomposition £(t) = 53j=i C;(£)® e jj we get the equivalent system 
ofSDEs 

dCiW = + f] C-(*-)d* + E ^Cj(*-)dW Q (t) 



t^2 71 



E ( -4= 4 - 1 ) o(*-)dJv«(*) + E E Ri k Ck(t-)dw*(t) 



a=d 1 + l v ^" ' o=l fc=l 



+ E E -7rr«i*c*(*-) -<*•(*-) djv^(t). (is) 

Q = d 2 + lfc=l \V^Q / 

As usual the solutions of SDEs with jumps are taken to be continuous from the 
right with left limits (cadlag processes); the notation i_ means the left limit. 

Remark 5. If some of the operators S in the jump part is zero, we eliminate its 
contribution by taking the corresponding Poisson process with zero intensity, so 
that it is almost surely for all times. In other words, if we have — for 
some k and some a > d,2, we take 4- 0- 

Assumption 3. We take a random normalized initial condition: £(0) = C° = 

n n 

EC°® e *> C° is ^-measurable, E® ||C°|| 2 ] = E E Q IK°|| 2 = 1- To repro- 

4=1 1=1 

duce the initial condition ([5} we ask also E Q []C°><C°l] = vW- Mean values of 
random operators are defined in weak sense. 

Equation (IT^I) is a particular case of the equations studied in Refs. |3QII40j . so, 
we refer to those papers for the properties of its solution, while all the results 
could be obtained by standard arguments in stochastic calculus and the Ito 
formula for continuous and jump processes summarized by the ltd table 

dW a (t)dW p (t) = S al3 dt, dW*(t)dW l p (t) = 5 af3 5 k idt, 
dN a {t)dN p (t) = 6 a pdN a (t), dNi(t)dN^(t) = 5 aP 5 l3 dN l a {t); ' 4 ' 
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all the other products are vanishing. 

Theorem 1 ( 30, Prop. 2.1, Theor. 2.4, Prop. 3.2]; gOJ Theor. 1.1, Theor. 
1.2]). Under Assumptions]]] and\3] the SDE (|12|) admits a unique (up to Q- 
equivalence) solution £(f), f > 0. Moreover, the mean state 
satisfies the master equation (O. 

n 

Finally, under the probability Q, the process p(t) :— ||C(f)|| 2 = IICiMI| 2 * s 
a non-negative (3^) -martingale with Q-mean 1 and if satisfies the Doleans SDE 

"" '/„(<) 



dKf)=K*-){E u «(*) dW '«(*)+ E (-^ 



- 1 (dN a (t) - X a dt 



7T12 n 



■EE«*)+ E E 

Q!=l fe=l Q = d 2 + lfc=l 



A* 



l] (diV Q fe (f)-A*df)], (15) 



w a (f) = 2Re(v>(f-) V a il>(t-j) = 2^fe(^(t_) i^j(t-) 



3=1 



3=1 



7 ;3 (f) = ||^(f-)|| 2 ^E|| L ^(*- 



3=1 



TTie process 



is defined by 



^(f)=ii^-)ir=E|<^(*-) 

3=1 



(16a) 
(16b) 
(16c) 
(16d) 

(17a) 

(17b) 



^ ( * ) = rer i/||cw|| ^°' 

V>, IIC(*)II =0, 

where ip £ CKs is a fixed vector of norm 1 / y^n. 

Remark 6 (A first unravelling). By the theorem above, Eq[|C(*))(C(*)|] satisfies 
the master equation © with initial condition 77(0) (Assumption [3]). So, we have 
77(f) = E Q [|C(f))(C(f)|], Vf > 0, and, by the discussion below Eq. (|T0]), we get 



m(t) =®®[Mt))(Ci(t)\], 



i = 1, . . . , n, t > 0, 



(18) 



which shows that £(f) is a pure-state unravelling of the solution of the Lindblad 
rate equation ©. 
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Remark 7 ( [30| Theor. 1.2]; [H Theor. 29.2]). The solution of the Doleans SDE 
(fI5]l is 



d 2 n 

+ EE 

mi ^ 

>< n 



/3=di + l 
m n 



a=l fe=l 
ft 



: (s)dW ra (s) 



v k a (s)dW*(s) 



2 Jo 



2 7o 



tiQ,(s) 2 ds 



^( S ) 2 ds 



(X -Ip(s))ds 



n 

re(0,t] 



Mi-l)AiV,(r) 



n n -p / (A*-w)d- n 

mi+l<=l ^ LJo J re(0,t] 



lW|(r 



where AN p {r,u) = N p (r,u)-N p {r-,u), Atf|(r,w) = N*(r,Gj)-N«(r-,oj). By 
the fact that a Poisson process has only a finite number of jumps in a compact 
interval, for every w only a finite number of factors contributes to the product 
over r in the representation above. 

Note that, if for some t, ui, j3, £ one has li(t,uj) — and ANi{t,uS) = 1, 
then p(T,u) = 0, VT > t. Similarly, Ip(t,u) = and ANp(t,u) = 1 imply 
p(T,w) = 0, VT>t. 



3.2 The generalized stochastic Schrodinger equation 

The final aim is to derive an equation for the normalized process (|17[) . This is 
based upon Ito stochastic calculus again and a Girsanov-type change of measure. 

Remark 8 (The change of probability measure). For for every T > 0, we define 
the physical probability P T over (fi,3V) by 

P T (A) =E Q [l A p(T))= I ||C(?»|| 2 Q(dw), VieJ r . (19) 

J A 

Note that P T depends also on £°, which we assume to be normalized in the 
sense of Assumption [3J The martingale property given in Theorem Q] ensures 
that the family of probabilities {P T , T > 0} is consistent, that is 

0<t<T, AeJt P T {A)=P t (A). (20) 

To obtain from (|2U)) the existence of a unique probability in the infinite 
horizon limit T — > +oo is a delicate problem and can be guaranteed only with 
respect to some sub-filtration composed by Borel standard cr-algebras [27J Sec- 
tion A. 5. 5]. 

It is important to note that the denominator ||C(t)|| in the definition of 
the processes ipkit) could indeed vanish as stated in Remark [7] But, by the 
construction in Remark[51 this happens with probability zero with respect to the 
new probability P T , while this is not guaranteed under the reference probability 

Q. 

The important consequences of this change of measure are the modification 
of the characteristics of the driving processes (t) and (t) (due to some 
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extension of the Girsanov theorem to the diffusive /jump case [H]) and the 
fact that ij)(t) satisfies a non linear SDE, the stochastic Schrddinger equation 

[22H2H1B]. 

Theorem 2 ( [3D1 Prop. 2.5, Thcor. 2.7]; [40, Prop. 1.1, Theor. 1.3]). Under the 
probability P T , </ie processes W a , W°, t <E [0, T], a — 1, . . . , d\, (3 — 1, . . . , d<x, 
k = 1, . . . , n, defined by 

W a {t) = W a (i) - / V a ( S )ds, = - / ^(S)ds, (21) 

Jo Jo 

are independent standard Wiener processes and the processes N a (t), Ng(t), t € 
[0, T], a = di + 1, . . . , mi, /3 = cfe + lj • • • , th-x, fe = 1, . . . , n, are counting 
processes of stochastic intensities I a (t) and Ig(t), respectively. 

Again under the probability ¥ T , the components of the process ip(t) satisfy 
in the time interval [0, T] the SDE 

di ^ 

dVi(*) = • ■ • , V»(*-))dt + E ( X « " 2 u «(«))^(i-)dW Q (t) 

a=l 

d 2 n 

+ E E ( R i k Mt-) - ^ «S(t)^(t_))dwj (t) 



d 2 n 

a=l fc=i 
mi 



+ £ tf^Hr "*('-)) dJV iW, (22a) 

a=d 2 +lfc=l V V 1 a\ t ) / 



1 _ TO1 

V j (Mt-),---,i>n(t-))=K j Tpj(t-)+2 E J a(*)V»i(t-) 

a=di + l 

ma ra 



mi n 1 / 1 \ 

2 E E^w-J^E^k-^aW 

Q= d 2 + lfc=l a=l ^ ' 

o E E < (*) - t «S(t)^-(t-) . (22b) 



d 2 n 

2 

a=lfe=l 

Note that the SDE (J^U) is non-linear in ip(t), because the quantities I a {t), 
Ia(t), v a (t), v a (t) are bilinear in tp(t) itself. Moreover, to consider (|22l) as a 
closed equation for tjj(t) poses interesting mathematical problems on the defini- 
tion of solution and on the meaning of uniqueness because the law of the driving 
noises depends on the solution tp(t) itself through the stochastic intensities 

la, It S3- 

Proposition 3 (A normalized unravelling). The solution of the Lindblad rate 
equation ([6]) can be expressed as the following mean with respect to the physical 
probability 

Vl {t) =E r r[\Mt))(Mt)\], i = l,-..,n, T>t>0. (23) 
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Proof. Let us introduce the set A t — {lu G O : cj)|| = 0}. Then, by the 
definitions of p(t) and ij)(t) given in TheoremQ] we have 

ICi(*)XC*(*)l = iA"|Ci(*)XCi(*)l = i Ai p(tM(t)){Mt)\ = P (t)\A(t)){A(t)\- 

By taking the Q-expectation and by taking into account Eq. ([TBI and the defi- 
nition of the new probability, we get 

Vi(t) = MKi(t))(Ci(t)\] =E Q [p(t)IV'i(*))(^(*)l] = %'[IV'i(*))^i(<)|]. 

Finally, by the consistency property (|2"0)) . we get (|23l) . □ 

This proposition gives an unravelling of the Lindblad rate equation ([6]) based 
on the components of the normalized vector ip(t). When d\ = m\ = and 
c?2 = 0, we recover the pure jump unravelling proposed in Ref. [34;. If the aim 
is only to simulate Eq. ([5]), a normalized pure state unravelling is much more 
efficient than a non- normalized one such as (|18[) [4J . The simulation techniques 
based on (|22]) with d\ = d^ = correspond to the Monte-Carlo wave function 
method started in Ref. |26) . while the case d\ = rti\ and e?2 = to 2 gives rise 
to simulations of diffusive type as in Refs. [48-50 . From the point of view of 
simulations, the fact that the starting point was Eq. ©, and not Eq. ([5]), has 
produced a more convenient unravelling with less noises (no dependence on the 
label j). 

4 Measurements and stochastic master equations 

In this section we face the problem of the measurement interpretation of the 
unravelling we have constructed. We introduce the notions of instruments and a 
posteriori states and we derive the non Markovian generalization of the stochas- 
tic master equations. 

4.1 Outputs and noises 

In the theory of measurements in continuous time [27, 30, 39 , 40 it is assumed 
that the output of the measurement is given by some components of the driving 
noises appearing in the SDE (Eq. (fl2l or (flU|) in our case); the law of the output 
in [0,T] is the physical probability (|1~9"|) . Not all the components of W and N 
have to contribute to the output. The role of some of the components of the 
noises could be only to perform the unravelling of some dissipative term. 

Let us examine first the components W k (t), a = l,...,d2, k = l,...,n. 
For t G [0,T], under the physical probability P T , from (J5T]) we get = 
Wq(*) + fo v a( s )ds', but, as one sees from Eq. (|16bp . v^(s) mixes different 
components of ip(t) and cannot be an observable, because it does not respect 
the superselection rule. In particular the mean value of W£(i) turns out to be 
E PT [W*(t)] = /„* E PT [«*(«)] ds with 

n 

E pT [v k a (t)] = E pt [v k a (t)} - 2^ReTr M {{W a k ® |e,)(e fe |) r,(t)} 

i=i 

and it involves the unphysical non-diagonal blocks Trc« {(1® l e i)( e fe|) Vi^)}- 
So, W£(t) cannot contribute to the output. 
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No problem of this kind arises for the other processes, as one sees from Eqs. 
(fTB) . The stochastic intensities I a (t), and the processes v a (t) do not mix 

different components of ij){t). However, if the counting process is detected 
we gain information on the block contributing to the emission (the block k), as 
one sees for instance from the mean intensity 

n 

^ [/£(*)] - E Tr « s { R i k * R i k vk(t)} . 

i=i 

If we assume that the index k is not physically observable coherently with the 
fact that the system state is the sum ((3]), the process is not observable by 
itself. However, there is no obstruction in considering as physically observable 
the counting process 

n 

M a (t):=5^JV*(t), a = d 2 + l,...,m 2 , (24) 
k=i 

whose stochastic intensity, under the physical probability, is Sfc=i-^aW- No 
problem arises on the observability of the other counting processes Np (f3 = 
d\ + l,...,mi), whose stochastic intensity under the physical probability is 
h(t). 

Let us stress that, under the reference probability Q, M a is a Poisson process 
of intensity 

n 

A « = X^a> a = d 2 + l,...,m 2 . (25) 
fc=i 

Let us consider finally the processes W a . At least in quantum optical sys- 
tems, observations with a "diffusive" character come out from heterodyne or 
homodyne detection and the involved operators must have an explicit time de- 
pendence due to the presence of the local oscillator [37J Chapt. 7]. We assume 
a very smooth time dependence which does not cause any essential change in 
the previous results. 

Assumption 4. For a — l,...,di, we assume the operators L 3 a to be time 
dependent and given by 

Ll(t) = hUtjU, , L{ e £(JC S ), \hi(t)\ = 1; 

the complex functions h J a (t) are continuous from the left. 

No time dependence is introduced into the master equations of Section [2] 
The explicit time dependence involves only the terms with dW a in Eqs. (1121) . 
(|T5| . (|22a[) and the third term in the right hand side of (|22b[) ; moreover, from 
Eq. (|16a[) . we get 



v a (t) =2Re^hi(t)(ip j {t-)\LiTp j (t-)). 
i=i 

The key result of the previous discussion is that, due to the mathematical 
structure and the meaning of the discrete label in the states, only the processes 
W a (a = 1, . . . , di), N/3 (j8 = di + 1, . . . , mi), M 7 (7 = d 2 + 1, . . . ,7712) can 
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be considered as possible components of the output. However, some of the 
components could represent pure noises, not observed quantities. So, we assume 
that only the first components are observed. 

Remark 9. Let us take d[ < di, m[ < mi, m' 2 < m^. We assume that the 
observed outputs are W a with 1 < a < d'i, Np with di + 1 < /3 < m\, M 1 
with c?2 + 1 < 7 < n^a • ^ some set of indices is empty, no component of 
the corresponding process is observed. The law of the output is the physical 
probability (fl9)) . Finally we denote by {St, t > 0} the augmented natural 
filtration generated by the set of the observed processes. 

4.2 The linear stochastic master equations and the instru- 
ments 

It is possible to have only the observed processes as driving noises in the dynam- 
ical equations, but for this we need to work with density matrices and trace-class 
operators. Let us introduce the positive trace-class operators 

cT(i):=E Q [|C(t)XC(i)||St], <T i (t)=E Q [\Ci{t)){Ci{t)\\9t\- (26) 

This means to take the mean on the non-observed components of the noises. 
Let us recall that £(0) is connected to the initial condition 77(0) (given in Eq. 
(fTU| ) by Assumption [3J By the fact that So is trivial we get 

(i(0) = r;(0) E S(5t), ^(0) = r?i(0). (27) 

Proposition 4. In the stochastic basis (f2, 5", St, Q), the operator valued process 
a(t) satisfies the linear stochastic master equation 

d[ 

da(t) = C[<r(t-)]db + J2 (V a (t)a(t-) + tr(t-)V a (t)*) dW a (t) 



a=l 



E ( K<7(t - )y " - a{t.) ) (djV a (t) - A a dt) 



+ E E " A ° -v(t-))(dM a (t)-A a dt), (28) 



where Ki{pi, . . . ,p n ) is defined by Eq. (|6b[) . L by ([8]), A Q by (|25|) . 

n 

V a (t) =Yy~Jt)t l a ® \e t )( ei \, a = 1, . . . ,di. 

i=l 

Given the initial condition, Eq. (|28p has pathwise unique solution. For the 
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components (the blocks on the diagonal) Eq. (|28[) reduces to 
da (t) = fCj(ai(t-), . . .,a n (tJ))dt 

d'l 

+ £ (hi(t) + hi(t)a 3 (t_)Li*) dW a (t) 

a=l 

+ E ( L> ^- )L « - (dN Q (t) X a dt) 



+ E (E R " ak{ l- )R * * - aj^)) (dM a (t) - A a dt) . (29) 

a=d 2 + l \k=l a I 

Proof. By applying the Ito formula to |C(*))(C(*)I an d, then, by taking the 
conditional expectation, we get the linear stochastic master equation for o~(t) as 
explained in [30l Sect. 4.2]. Existence and uniqueness of the solution of Eq. ((28)) 
is given in [301 Prop. 3.4]. Equation (|29p is obtained by direct computations. □ 

Remark 10. Let us consider now the physical probability introduced in Remark 
[51 The probability density of the restriction of P* to St with respect to the 
reference measure Q on Si is 

n 

p s (t) = E Q \p(t)\9t] = TxMt)} = ElW^(*)}- (30) 

The density pg is a S-martingale under Q and the restrictions of the physical 
probabilities are consistent. 

In the axiomatic formulation of a quantum theory, measurements are rep- 
resented by instruments, which give the probabilities and the states after the 
measurement (o posteriori states). As in Refs. [301140] . we put 

l t (F)[r l (0)}=E Q [l F <j(t)}, F€9t, i7(0)eS(5£). (31) 

By linearity we extend 1t{F) to the whole T(!K) and we get an instrument with 
value space (f2, St), which means that It(F) is a CP map from T(J{) into itself 
for all F € St, it is a strongly er-additive measure as a function of F and It(fi) 
is trace-preserving. 

Remark 11. Let us particularize the definition of instrument in the enlarged 
space to our case. We define 

T t (F)[r )1 (0),...,r ]n (0)}=E®[l F a i (t)}, F e St, (32) 

for all 77(0) satisfying the superselection rules. With the notations of Remarks 
ElandHlwe have (2j(F), . • • ^t( F )) = V oI *( F )| e ( x . T(Ms) ) ■ This is an instru - 

ment with the same value space as before, but made up of maps on C(X; T(J{g)) . 
Finally, by defining 

n 

lf(F) = E^'(n (33) 

3=1 
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we get an instrument with value space (fl,St) made up of CP maps from 
C(X;T(!Ks)) into T(!Ks). Moreover, the connection with the various dynamical 
maps introduced in Remark 2] is given by I{(0) = T(t), 

(^(^),-,W))=^°T(i)| e(x;TC;Ks)) , ^(P)=P 8 oPoT(t)\ e{xinxB) y 

The instruments give the physical probabilities once one has the pre-mea- 
surement state. In our case we have, V-F S St, 

n 

V\F) = Tr x {It(F)[ri(0)]} = £ TW^JMO), . . .]} 

i=i (34) 

= Tr Ms {li(F)[r h (0),...}}. 

This equation says that the probabilities P*, introduced before by starting from 
some stochastic differential equation, can be obtained also from instruments; so, 
the axiomatic structure of a quantum theory is respected and the interpretation 
as physical probabilities is justified. 



4.3 The a posteriori states and the stochastic master equa- 
tion 

The instruments give also the a posteriori states, the conditional states after the 
measurement. Let us recall the definition in the case of It; in the other cases 
the definition is analogue. The a posteriori state for the instrument It and the 
pre-measurement state 77(0) is the §(5f)-valued random variable p(t) such that 

X t {F)[r l (0)]=E rt [l F p(t)}, VFeS t . 

By taking into account that the density of P 4 with respect to Q is the trace of 
a(t) and how I t (F) [77(C))] is defined in terms of cr{t), we get easily 

p(t) = T jff W} =MM*))M*)IN- 

The components of p(t), which are 

IQJ = 

Tr M {a(t)} 



Pi(t)=Er[\Mt))Mt)\\9 t ] = °fL x , i = l, 



give the a posteriori states for X\: 

2i(F)[r/ 1 (0),..., ?? „(0)]=E pt [l F p i (t)]. (35) 

Note that we have also Ep^/o^t)] = r/i(t). Finally, by taking the sum over i in 
Eq. (|35[) . we get the a posteriori states ps(t) = Pi(t) for Xf . 

On the other side, the states r/(t), rji(t) are called the a priori states, due to 
the fact that these states are the averages of the a posteriori states and that 
they are interpreted as the states to be assigned to the system at time t when 
the result of the observation is not taken into account. 
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Remark 12 (The stochastic master equation). For a — 1,...,^, /3 = d% + 
1 , . . . , m' 1 , 7 = d,2 + 1 , . . . , m' 2 , let us define 

n 

m a {t) =: Ept [» a (t)|St] =2Re^/ii(i)Trjc s |L*p rf (t_)} , (36) 

i=i 

n 

J^(t) := Ep. [Ip(t)\9t] = E Tr ^ > ( 37 ) 

3=1 

n n 

J 2 (i) := ^E pt [/ 7 fe W|St] = £ {<*<p fe (^)} • (38) 

fc=i j,fc=i 

Then, by stochastic calculus, under the new probability P T and for £ € [0, T], 
we get the equation for p(t) [301 Rem. 3.6] and, then, the stochastic master 
equation for the components 

dp j (t) =JCj(pi(t-),...,p n (t-))dt 

+ £ (V a (i) + hi(t) Pj (t-)Li* - m a (t) Pj -(t_)) dW" a (t) 

a=l 

+ E " -ft(t-) (d^( f )-4(t)dt) 

+ E (E ^tl^ - Pi(*-)1 (d^rW ~ ■ 09) 

7 = d2 +i \fe=i "M^ / 

The processes lV a are independent standard Wiener processes, iVg(t) is a count- 
ing process of stochastic intensity Jl(t) and M 7 (t) is a counting process of 
stochastic intensity 

5 A two-level system in a structured bath 

To give a simple, but concrete example of the theory we have developed and to 
have a first idea of the effects on physically measurable quantities, here we study 
a model of a two level atom in contact with a non-trivial structured reservoir and 
we compute the heterodyne spectrum of its emitted light. This is a modification 
of a model |15 | I16 | I3"3] which could represent the dynamics of a single qubit in a 
non Markovian environment or the dynamics of an optically active molecule, as 
the fluorophore system, in a local nano-environment [2~i] . 

We consider a two-level system in contact with a two-band reservoir; so, 
!Ks = C 2 and n — 2. Let a z , a± be the usual Pauli matrices; then, P + = a + a- 
is the projection on the excited state (J) and P_ = cr_<7 + the projection on 
the ground state (°) . Here we give the mathematical model, while the physical 
interpretation is given when we write down the various dynamical equations. 
By using the notations introduced in Assumption [2] and Section 13.11 the model 
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we consider is denned by the following choices: 

d 1 =mi = 2, d 2 =0, TO 2 = 2; H* = y cr z , u; 4 > 0, i = 1, 2; 
i?f=0, flf = V7r^-, i?l 2 = V7Ja+, 7l >0, 7 = 1,2, 
i$=0, i?2 2 = 0, i?^ = V7o^l, >^>0, 7o >0; < e < 1, 

L\{t) = L\ (t) = c M ^a- , L\ (t) = L\ (t) = e ivt ^(H <r_ , ^1. 

(40) 

The driving processes in the linear SDEs are the standard Wiener processes Wi , 
W2 and the Poisson processes iV-J, Nf, with intensities A}, \\ , X\\ all these 
processes are independent. According to Remark O take \\ \. 0, so that TVf is 
almost surely and we can set diVf (t) = 0. 

5.1 The Lindblad rate equation and the equilibrium state 

First of all let us write down in the concrete case introduced above the Lindblad 
rate equation ^ 

^ m (t)=lC 1 ( m (t),i l2 (t))= 1 o (a-r h {t)a + -h L P + ,r n {t)} 

~.| i^x 
+ I2<r+V2{t)v- - y{-P+>»7i(*)} -70^1 (*) - — [o-2,f?i(*)] » (41a) 

^i»(t) =AC 2 (»7i (*),%(*)) =7o fff_»ft(t)o- + - i{P+,7 /2 (t)} 

+ 7l( T_ry 1 (t) ( T + - ^{P_, 772(*)} +7o»7i(t) - ^ Wz,m(t)} ■ (41b) 

The model of Refs. [T^ll^Tj corresponds to 70 = 0, x = 0, u>i = o; 2 ; moreover, 
the rotating framework is used, so that the terms with u!i disappear. In Refs. 
[THUD t ne case w i 7^ w 2 is allowed and it is explained by different energy 
shifts induced by the two bands of the environment. So, we have a two level 
molecule with two resonance frequencies due to the structured environment. 
The terms with 71 and 72 represent the molecular transitions induced by the 
environment and concomitant with transitions between the two bands of the 
nano environment. 

Reference [23] studies the stimulated fluorescence light under laser excitation 
of the molecule; the treatment is based on the quantum regression formula. 
Instead, our aim is to study the spontaneously emitted light and to this end 
we have added the first term in both equations, the one with 70, which is an 
explicit spontaneous emission term. 

To have emission without stimulation by external light, we need some thermal- 
like excitation. To get this effect we have added the terms with jo>c. This is 
the simplest modification giving rise to a non trivial equilibrium state. 

If we write Eqs. (|4ip in terms of the matrix elements of 771 and 772 we get two 
decoupled systems of equations: the system for the coherences (the off diagonal 
terms) and the system for the populations (the diagonal terms). Firstly, one 
checks easily that the coherences decay exponentially to zero. On the other 
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side, the system of equations for the diagonal terms turns out to be equivalent 
to a 4-state, irreducible classical Markov chain. If we denote by 1 + , l - , 2 + , 2~ 
the four states, the transition rates different from zero are 70 for the transition 
1+ -> 1~, 71 for the transition 1+ -> 2~, 70^ for 1+ -> 2+, 72 for 2~ -> 1 + , 70 
for 2 + — > 2~, 70 x for l - — > 2~. From the graph of this finite-state Markov chain 
we see that it is irreducible; then, there is a unique equilibrium distribution, 
computed below, and it is a global attractor. 

Equilibrium state. 

The Lindblad rate equations (141[) admit a unique equilibrium state 7]i(oo) = 
lim i _ ) . +00 7]i(t), i — 1,2, which can be easily computed. It turns out to be given 

by 



v%(oo) = Pi {4 p + + z i p -) > z i ■= 1 - z ii 



z := 



1 + Xi ' 



72^ 

■ 7— — r, Pi:=P, p 2 :=l-p, 

7i +7o(l + 

72(1 + x) 



72 + x (7o + 72 + 7i) + x 2 (70 + 72) 



Let us note that we have (1 — p)z^ = xpz^ . By recalling that the system state 
is the sum of the components ©, we get that the average equilibrium state of 
the two-level system is 

% q = + ^2(00) = pxP + + (1 - px)P_ . 

5.2 The stochastic Schrodinger equations 

The 1SSE (JI3J corresponding to the choices (gOJ) is 



dCi(t) - (^ + ^) Ci(t-)dt-Ci(t){dNl(t)+dN^t)) 



[,jf 2 a + C 2 (t-)-Ci(t-))dNf(t) 



■-2 



+ C a (*_)dt + U^Cx{t-) - Ca(*-)) dJVaHt) 
C 2 WdA r 1 2 W + ff -Ci(M - Ca(*-)) dJVj(t) 



dC 2 (t) = x 



+ e il/t V7^^-C 2 (i-) (VidWi(t) + vT^dW 2 (i)) , 

where A = A} + \\ + A^ and 

^ 1 = _^i 70+71 p 7ox K 2 = _^2 To 72^ 

2 2 2 2 2 2 

Note that the Wiener processes W% and W2 appear always in the combination 
y/eWx(t) + \/l — e W 2 {t), which is again a one-dimensional standard Wiener 
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process. The reason for the introduction of two components is that the diffusive 
term represents the emitted light, which we have divided in two channels: chan- 
nel 1, represented by W\, contains the light reaching the heterodyne detector 
and channel 2, represented by W2, contains the lost light. The proportion of 
lost light is 1 — e. 

Finally, by Eqs. (j2"2")l . the SSE for the normalized vectors is, under the phys- 
ical probability, 

<tyi(i) = yi(V>i(t_),V2(t-))dt- ViWCdiV^t) +dN'(t)) 



-ih(t-))dNt(t) 



h+Mt-)\\ 

" - ' ' ' ' ' " ' Mt-)) (t) -Mt- )dN? (t) 



.\W-Mt-)\\ 

+ V70 U vt a-Mt-)- ^«(t)^(*-)J (VedW^t) + VT~edW 2 (t)^j 
where W is the new Wiener process introduced in (j2"Tj) . and 



V\{t) = \^m v (t), v 2 {t) = V7o(l - e)v(t), 

2 

;(*) = 2^Re(e i ^ fe (*-)|<7-^(t-))), It(t) = x l0 Wi(t-)\\ 
fc=i 

/J(t) = 71 \W-Mt-)\\ 2 , /i 2 W = 72 h+Mt-)\\ 2 , 



5.3 The stochastic master equation and the heterodyne 
spectrum 

In the situation we are considering the band transitions cannot be monitored. 
We take under observation the system by collecting part of the emitted light 
in an apparatus performing heterodyne detection. In this detection scheme the 
received light is made to interfere with some monochromatic light of frequency 
v; to a certain extent, this frequency can be varied. As we have said, it is 
W\ which represents the light reaching the detector; moreover, the (stochastic) 
output J{t) of the detector is some smoothed version of W\ [371 Chapt. 7], say 

J{t) = Vfc f e-^-^/ 2 dWi(s), k > 0. (42) 
Jo 
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To take into account that only W\ is observed we use the notation of Remark 
[5] and we take d[ — 1, = d\, m' 2 = 0; recall that we have d\ — mi = 2, 
c?2 = 0, 77i2 = 2. Then, all the sums with jump processes disappear from the 
stochastic master equations (|29p and (|39l) . The linear stochastic master equation 
((29)) becomes 

doj-(t) = ^((7i(i),o- a (t))df + (e'^-o-^t) + e^V^i)^) dWi(i), 

where the JCj are the operators appearing in the Lindblad rate equations PTj). 
The corresponding non linear stochastic master equation (|39[) for the a posteriori 
states pj(t) — <Jj(t)/ r Iv'K s {ai{t) + cr 2 (t)} turns out to be 

+ (e iut a- Pj {t) + e-' wt P] {t)a+ - m(i) Pj {t)) dWi(t), 

m(t) = 2 Re (e 1 "' Tr Ms {a_ Pj (t)}) . 

The power of the output current produced by the detector is proportional 
to J(t) 2 and the mean power at large times is proportional to 

P(v) = lim E pt [J{t) 2 }. (43) 

t— > + oo 

The limit is to be taken in the sense of the distributions in v. 
By using (|42l) we get 

I e fcs/2 d^i(s) / e fcr / 2 dWi(r) 



Ept[J{t) J ] = ke- k % 

So, to obtain an explicit expression for the power first of all we need to compute 
the second moments of the Wiener type integrals L e ks / 2 dWi(s) under the 
physical probability. The autocorrelation function of W\, from which such a 
moments follow, can be obtained by differentiation of the so called characteristic 
operator (the Fourier transform of the instruments) |27l Proposition 4.16]. The 
formula valid for the Markov case needs only to be expressed by using the 
diagonal blocks; from [27J Eq. (4.47)] we get 

E pt [J(t) 2 ]=k /V fe ('- S ) ds + 2fc 70 e fds [ S dre- k ^- s ^ 2 e~ k ^/ 2 
Jo Jo Jo 

2 

x £ Tr Ms {(e-V_+e- i "V + )^( S -r) [e ws a. m {r) + e'^r)^] } . 

(44) 

By 2i=i J~ij (*) i T j] i * = 1) 2) we denote the solution of the Lindblad rate equation 
(|4ip with initial condition (n, r 2 ) at time 0. Then, the computations needed to 
obtain P(v) are long, but similar to the ones in [27J Sect. 9.1]; we give only the 
final results: 

P(v) = 1 + 4%eE(v), (45) 
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E(V) = 2 7o 



2 7 <r 



(1 - p) z^T 2 - pzf w [r 2 (72 - 71 - 270^) + 4 (w 2 - wi) {v - lo 2 )\ 



pz\ {[1 + W (72 - 71 


- 2 7o >ir)] Ti + 4w 


(w 2 - wi) - Wi)} 


7T 


4(^-w 1 ) 2 +r i 2 





(46) 



where p, zj~ have already be defined and 

Ti := 70 + 7i + 2 7o>< + fc, T 2 := 70 
270^ 



72 



w :- 



4 (wi — oj 2 ) 2 + (72 — 7i — 27ox-) 2 



In Eq. (|4"5|) the term 1 is interpreted as the shot noise due to the local oscillator 
and as the heterodyne spectrum. Note that the widths Tj contain some 
dynamical parameters and the instrumental width k. 

By its definition, we have P{v) > 0, while the positivity of the spectrum 
E(i/) is not obvious. However, one can check that it is possible to rewrite T,(y) 
in a form from which its positivity is apparent: 



S(z/) = Z?7qX'< 



7o(l + x) +71 + k 
A{u- ujrf + 1\ 2 



x (72 + k) 



4(1/ — W2) - 
j x(T 1 



r 2 2 



D = 



i(v-u 1 y + r^J [4 {v - io 2 y + r 2 2 

2/tt 



(47) 



1 + ^71/72 + x(l + x)(l +70/72)' 

Note that the heterodyne spectrum, for spontaneous emission in our model, 
contains information on the dynamics: all the dynamical parameters, due to 
the structured reservoir, determine the form of the spectrum. In particular, we 
have a double peaked structure only if wi and w 2 are sufficiently different and 
this difference can be only due to the band structure of the bath. 
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